Pinned states in Josephson arrays: A general stability theorem 
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Using the lumped circuit equations, we derive a sta- 
bility criterion for superconducting pinned states in two- 
dimensional arrays of Josephson junctions. The analysis ne- 
glects quantum, thermal, and inductive effects, but allows 
disordered junctions, arbitrary network connectivity, and ar- 
bitrary spatial patterns of applied magnetic flux and DC cur- 
rent injection. We prove that a pinned state is linearly stable 
if and only if its corresponding stiffness matrix is positive 
definite. This algebraic condition can be used to predict the 
critical current and frustration at which depinning occurs. 

PACS Numbers: 74.50. +r, 05.45.+b, 03.20.+i, 85.25.Cp 



Collective pinning occurs in a wide variety of coupled 
physical systems. Examples include vortices in Type-II 
superconductors, cracks and dislocations in solids, and 
charge-density waves in quasi-one-dimensional metals.EJ 
In each case, when the system is subjected to an exter- 
nal constant drive, it remains motionless until the drive 
exceeds a critical value (the depinning threshold) after 
which the system begins to move. The pinning is collec- 
tive in the sense that it involves interactions among many 
coupled subsystems, typically in the presence of disorder. 
Hence it is often difficult to predict the depinning thresh- 
old theoretically. 

Here we study collective pinning for a relatively 
tractable class of model systems: two-dimensional (2D) 
arrays of Josephson junctions. Besides their technologi- 
cal applicationsp Josephson arrays can be used to explore 
fundamental questions in statistical mechanics (such as 
phase transitions), and in nonlinear dynamics (such 
synchronization and spatiotemporal pattern formation) 
In addition, they have been proposed as clean models 
for layered and granular high-T c super conductors. As 
such, their depinning could be relevant to the under- 
standing of the onset of resistance in the current-voltage 
characteristics of high-T c samplcs.13 

Sfjyexal advances hasie occurred recently in the numeri- 
calO'lltl and analyticaO~E3 investigation of 2D Josephson 
arrays, thanks in part to an influx of ideas from nonlin- 
ear dynamics. In this paper, we analyze depinning in 2D 
arrays from this perspective. Using a compact matrix 
notation, we show that the linear stability problem for 
pinned states can be mapped onto the classical mechani- 
cal problem of small oscillations in a network of coupled, 
damped linear oscillators. The results apply to 2D ar- 
rays of any given topology. There are also no restrictions 
on the capacitances, resistances, or critical currents of 



the junctions, nor on the spatial patterns of DC current 
injection and applied magnetic flux. Our main result is 
that a pinned state is stable if and only if its correspond- 
ing stiffness matrix K is positive definite. This matrix K 
changes with the pinned configuration and depends on 
the connectivity and disorder of the array. A corollary is 
that any pinned state with all phases \<pi\ < n/2 is guar- 
anteed to be stable. We also prove that depinning can 
never occur via a Hopf bifurcation; only zero-eigenvalue 
bifurcations are possible. 

Our analysis is based on several simplifying assump- 
tions. First, we neglect thermal fluctuations; that is, we 
assume zero temperature. Second, we assume that the 
superconducting islands in the array are large enough 
that quantum (charging) effects are negligible. Thus, 
the phase 9i of the complex macroscopic wavefunction 
at each island is a well-defined classical variable. Third, 
we assume that the junctions between islands are small 
enough that they can be approximated as lumped ele- 
ments. Therefore, the junction between two islands I 
and to can be described by a point gauge-invariant phase 
difference 
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where A is the total magnetic vector potential and 
$o = h/(2e) is the quantum of magnetic flux. Fourth, 
we modeLeach junction by the standard RCS J equivalent 
circuitcltj with superconducting, resistive, and capacitive 
channels in parallel. Then the junction dynamics obeys 
a damped driven pendulum equation 
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with effective mass //$ = <I>oCi/(27r/ c o), damping 7, = 
$o/(27T-R,/ c o), and restoring strength <r)i = I cl /I cQ . The 
capacitance Ci, resistance and critical current I C i 
are fabrication- and material-dependent parameters that 
characterize junction i. The drive is given by the nor- 
malized current i\, measured in units of I c q = (Id), the 
average critical current of the junctions in the array. 

In dealing with. arrays, it is useful to introduce a vector- 
matrix notationjj'otilo where the variables are now vec- 
tors of three types: node vectors of dimension n (e.g. 
8); edge vectors of dimension e (e.g. <f> and i b ); and cell 
vectors of dimension c defined at each plaquette. (More 
precisely, n is the number of independent-nodes, after one 
node is grounded and taken as reference. o) The edge and 
node variables are related through an e x n edge-node 
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connectivity matrix A that encodes the topology of the 
array, including its boundary conditions such as the pres- 
ence (or absence) of edges. Similarly, an e x c edge-cell 
matrix B transforms between edge and cell variables, in 
what amounts to taking a discrete curl. 

Within this framework, the nonlinear constitutive 
law (|J) can be compactly written as 

fi(f> + 70 + r\ sin 4> = i h (3) 

where fj, = diag(/ij), and 7 and rj are similarly defined 
diagonal matrices. Each junction is allowed to have a 
different capacitance, resistance, and critical current, as 
recorded in the matrices 7, and 77. 

When junctions are interconnected to form a network, 
there exist topological constraints which can be expressed 
in terms of the connectivity matrices A and B~-. First, the 
currents must satisfy Kirchhoff's current lawii-3 

A T i h = i cxt , (4) 

where the vector i cxt gives the balance of normalized 
current at each node, and reflects the particular scheme 
of current injection/extraction for each experimental de- 
vice. For instance, in the usual experimental setup, 
where a uniform DC current I^c is injected (extracted) 
at the bottom (top) nodes, all the components of i ext 
will be zero except those at the bottom (top) boundary, 
which will be equal to idc ( — idc)- Our analysis, however, 
is valid for an arbitrary injection scheme, as long as the 
bias currents are time-independent. 

The second topological constraint is the flux quantiza- 
tion in each cell of the array. We assume the simplest 
case where all self-fields due to inductance effects are ne- 
glected. Then the flux quantization is given by 

B T ct) + 2ttF = 2n( = 0, (5) 

where £ is a cell vector of integers (topological vorticities) 
that have no dynamical relevapee, and can be redefined as 
zero with no loss of generality.liil The cell vector F records 
the external flux through each plaquette, measured in 
units of the flux quantum. In experiments, the external 
magnetic field is often spatially uniform across the array. 
Then F is a constant vector with value / = < £ , cxt/ < l > o- Our 
analysis holds more generally for any time-independent 
spatial pattern of applied flux. 

For the no-inductance case assumed here, the transfor- 
mation (Q) between junction and island phases is given 
in vector form by 

(j> = AO - ip (6) 

where <p is a time-independent edge flux vector, fixed by 
our choice of gauge but subject to B T ip = 2irF, which 
follows directly from (||) noting that BJ.A = 0, from 
the definition of the topological matriceso From (g), (jij) 
and ([}]), we obtain the governing vector equation of the 
system: 



A T fiA 9 + A T ~/A 9 + A T j] sm(A6 - <p) = i cxt . (7) 

From now on, we focus on the pinned states of the 
array. These correspond to static configurations 9* of (|?|), 
given implicitly by 

A T rj sm(A8* — ip) = z cxt . (8) 

Typically this nonlinear algebraic system (|^) has multiple 
solutions. Each solution depends parametrically on the 
external current vector i oxt and the applied flux vector 
F. (In the usual experimental setup, these are deter- 
mined by the scalars idc and /, respectively.) As i cxt or 
F are varied, the linear stability of a given static config- 
uration 9* can change. This signals the transition to an- 
other state of the system. If the new state is still pinned, 
the transition corresponds to a static rearrangement of 
phases and currents; on the other hand, if the new state 
is time-dependent, it corresponds to depinning and the 
onset of resistance. (Because our analysis is local, it can- 
not distinguish between these two types of transitions.) 

To study the stability of the pinned states, let = 
9* + a where a is a small perturbation. Linearizing (Q) 
about 9* yields 

Ma + Ga + Ka = 0, (9) 

where 

M = A T nA, G — A T jA 1 K = A T r]C*A (10) 

are the mass, damping, and stiffness matrices, respec- 
tively, and 

C* = diag(cos0*) (11) 

is a diagonal matrix of the cosines of the phases of the 
given static configuration. Both M and G are symmetric, 
positive definite matrices, since A is a topology matrix 
and fx and 7 are diagonal matrices with positive masses 
and damping coefficients on the diagonalo However, K 
is not necessarily positive definite since the cosines on the 
diagonal of C* are not necessarily positive. We stress 
that the stiffness matrix K is different for each pinned 
state, and it also changes parametrically with the exter- 
nally tunable parameters. 

Equation (j9j) is familiar from the classical mechanical 
problem of small oscillations-jn a network of coupled, 
damped harmonic oscillatorsJla But the present stability 
problem is not as trivial as it might seem. Ordinarily one 
assumes that K is positive definite, but that need not be 
true here. Also, recall that when damping is present, 
normal modes cannot be used to decouple the system; 
in mathematical terms, one cannot simultaneously di- 
agonalize the three symmetric matrices M, G, and K. 
Therefore we analyze (|^) from first principles. 

A given pinned state is linearly stable if and only if the 
perturbation a(t) decays to zero for all initial conditions. 
Equivalently, all the eigenvalues of ([)]) must have strictly 
negative real parts. The characteristic equation 



2 



det(A 2 M + XG + K) = 



(12) 



cannot be solved explicitly, but one can still extract 
useful information about the eigenvalues, as follows. 
Suppose that ( |l2| ) holds for some A. Then there ex- 
ists a (possibly complex) eigenvector x ^ such that 
X 2 Mx + XGx + Kx = 0. Multiplying on the left by the 
complex conjugate transpose yields 



X 2 m + Xg + k = 0, 

where m = x'Mx, g = x'Gx, and k 
that depend on x. Thus, 



A 



-g ± \J g 2 — 4km 
2m ' 



(13) 

x' Kx are scalars 
(14) 



The key point is that m > and g > for all x, since M 
and G are real and symmetric (hence Hermitian) positive 
definite matrices. On the other hand, K is not necessar- 
ily positive definite, so k can have either sign. If k > 0, 
there are two subcases: if g 2 — 4km < 0, the eigenval- 
ues are complex conjugates with Re(A) = —g/(2m) < 0; 
otherwise the eigenvalues are both real and negative. In 
either case, the eigenvalues for k > lie in the left half 
plane and therefore correspond to stable modes. On the 
other hand, if k < 0, then A_ < 0, A + > so the A+ 
mode is unstable. Finally, if k — 0, then A- < 0, A+ = 0, 
and the A+ mode is neutral. 

An important qualitative conclusion from these formu- 
las is that any eigenvalue of ( |l2| ) must be either pure real, 
or complex with strictly negative real part. In particular, 
pure imaginary eigenvalues are forbidden. An immedi- 
ate consequence is that pinned states can never undergo 
Hopf bifurcations; depinnme can occur only through 
zero-eigenvalue bifurcationsc3 such as saddle-node, tran- 
scritical, and pitchfork bifurcations. 

We now prove the main result: a pinned state is lin- 
early stable if and only if K is positive definite. To prove 
the "if" direction, suppose that K is positive definite. 
Then k > for all eigenvectors x. From Eq. ( |l4| ) above, 
Re(A) < for all A and, hence, the pinned state is linearly 
stable. 

To prove the "only if" direction, it is equivalent to 
prove its contrapositive, i.e., we assume that K is not 
positive definite and show that the pinned state is not 
linearly stable. There are two cases. If det(K) = 0, then 
A = is a solution of (|l2|), by inspection. But A = 
corresponds to a neutral mode, not a decaying mode as 
required for linear stability. Next suppose det(if) ^ 0. 
We outline a homotopy argument which proves that Jl^ ) 
has a root A > 0. The strategy is to start with the un- 
damped problem, where it is easy to show that there is an 
unstable mode if K is not positive definite. Then we con- 
tinuously deform the undamped problem into Eq. (12), 



and show that the unstable eigenvalue remains unstable 
throughout the deformation. More precisely, consider the 
one-parameter family of equations 



where < p < 1 is a homotopy parameter. At p — 
0, Eq. (|l5|) corresponds to an undamped system, and 
normal modes can be used to show explicitly that ( |l5| ) 
has an eigenvalue A(0) > 0. As p varies continuously from 
to 1, this eigenvalue traces out a continuous curve X(p) 
in the complex plane. The curve starts on the positive 
real axis since A(0) > 0, and it must stay there for all 
p because any eigenvalue in the right half plane must be 
pure real, as shown by (|l4"|). Moreover, the curve cannot 
cross through the origin; from (p"5|), X(p) = for some p 
would imply dct(if) = 0, contrary to assumption. Thus 
X(p) > for all p. Setting p — 1 yields the desired result 
that (12) has a root A > 0. 



One consequence of this theorem is an implicit formula 
for the stability threshold of a pinned state 8*. As we 
vary the applied current or magnetic field, 9* and its as- 
sociated matrix K will change. The theorem implies that 
9* loses stability precisely when K — A T r\C* A ceases to 
be positive definite. This threshold is reached when the 
following algebraic condition is satisfied for the first time: 



det(if) = det(A T r]C*A) = 0. 



(16) 



Hence the stability threshold for 9* is determined exclu- 
sively by the array topology, by the injection scheme and 
bias current (through i cxt ), by the applied magnetic field 
F, and by the disorder in the junctions' critical currents 
(via the matrix rj). On the other hand, it does not de- 
pend on the mass (capacitance) and damping matrices M 
and G. This means that overdamped and underdamped 
systems have identical depinning thresholds. 
Another corollary is that if 



cos < 



>0, Vt 



(17) 



det(X 2 M +pXG + K) = 



(15) 



then that configuration is stable. This follows from the 
fact that the diagonal matrix r\C* of such a configuration 
is positive definite; therefore K is also positive definite.Ed 
On the other hand, since K can be positive definite even 
if C* is not, ( [l7|) is only a sufficient (but not necessary) 
condition for the stability of a pinned state. 

The constraint ( |l7| ) has a clear physical meaning for 
a single, isolated Josephson junction. Recall that as the 
bias current is increased from zero, a single junction re- 
mains pinned un±il_j^) = n/2, at which point it depins to 
a running mode.HHj Extrapolating naively from a single 
junction to an array, it is tempting to conjecture that 
an array should depin when its "most unstable" junction 
first reaches <j> — ir/2. Note, however, that this heuristic 
depinning criterion is equivalent to det(C*) = 0, rather 
than the rigorous condition det(K) = 0; therefore, it is 
not exact. Nevertheless, for the specific case of a ladder 
array with square plaquettes and perpendicular current 
injection, we have shown elsewhere£L3 that it can provide 
a good approximation to the true depinning threshold. 

The algebraic condition ( |l6| ) can be used to ease the 
numerical determination of the depinning threshold for 
2D arrays. For instance, the depinning current is usually 
obtainedQ through dynamical simulations that resemble 
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the actual experiment: the current is ramped up adiabat- 
ically and the circuit differential equations are numeri- 
cally integrated until a running solution appears. In con- 
trast, we solve (|8|) and (|l6| ) simultaneously to determine 
the critical current and the bifurcating phase configura- 
tion as functions of all the other parameters. This purely 
algebraic calculation can be done by Newton's method 
or some other rootfinding scheme. The results coincide 
with those found dynamicallyJij 

Another theoretical approach to depinning uses ther- 
modynamic and quasistatic calculations of pinned 
states.tZrt3 One can show that the condition ( filf ) is 
strictly equivalent to finding the point at which a given 
static configuration ceases to be a minimum of the po- 
tential energy 

V = -0U ext -Tr(r)C). (18) 

Thus a soft-mode conditional rigorously predicts de- 
pinning, while the criterion based on maximizing the 
quasistatic current induced by twisted boundary condi- 
tionsEa is only approximate.!] Note also that, although 
stable static configurations correspond to local minima 
of V, we do not attempt here to obtain the absolute 
minimum of the potential energy. This problem would 
require global optimization methods, such as simulated 
annealing. 

Our results open several promising lines of research. 
First, our analytical framework facilitates exploration of 
the effects of network connectivity on the depinning of 
Josephson arrays. The implicit condition (jl(^) can be 
turned into explicit, testable predictions of the applied 
current and frustration at which depinning should oc- 
cur. It may be possible to obtain analytical results for 
square and triangular arrays of identical junctions, pQ*-i 
haps along the lines of recent work on ladder arrays .113 
Second, one should also try to take self-fields into ac- 
count. Preliminary results suggest that the formulation 
given—here can be generalized to include inductance ef- 
fcctsJIS Finally, it is important to study more quantita- 
tively how disorder affects the stability of pinned states, 
both as the inevitable result of fabrication irregularities 
and as a design tool to manipulate the response of the 
network in a controlled fashion. 
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